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^ Abstract. Let Z be a positive integer and s = (si, . . . , si) be a sequence of positive integers. In this 

paper we shall study the arithmetic properties of multiple harmonic sum H(s;n) which is the n-th 
partial sum of multiple zeta value series C( s )- We conjecture that for every s and every prime p there 
are only finitely many p-integral partial sums H(s; n). This generalizes a conjecture of Eswarathasan 
and Levine and Boyd for harmonic series. We provide a lot of evidence for this general conjecture as 
well as some heuristic argument to support it. This paper is a sequel to Wolstenholme Type Theorem 
. for multiple harmonic sums, Intl. J. of Num. Thy. 4(1) (2008) 73-106. 

Resume. Soit I un entier et s = (si, . . . ,si) une sequence d'entiers positifs. Dans ce document, nous 
etudierons les proprietes arithmetique de sommes harmoniques multiples H(s;n), qui est le n-me 
somme partielle de la valeur de la serie multiple zeta C( s )- On conjecture que pour tout s et de tous 
les premiers p, il n'y a que de nombreux finitely p-partie integrante sommes H(s,n). Ceci generalise 
une conjecture de Eswarathasan et Levine et Boyd pour la serie harmonique. Nous fournissons beau- 
coup d'elements de preuve pour cette conjecture generate ainsi que certaines heuristiques argument 
\Q • soutenir. Ce document fait suite a Wolstenholme Type Theorem for multiple harmonic sums, Intl. J. 

of Num. Thy. 4(1) (2008) 73-106. 
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O ' 1 Introduction 



In recent years there is a revival of interest in the multiple zeta values defined by 



C(s) :=t(s 1 ,...,s l )= J2 k i Sl --- k 

0<ki<---<ki 



—si 
I 



for s — (si, . . . , si) £ N l , where N is the set of positive integers (see, for e.g., [IIJ [13l [19] ) . In physics, 



X. 

however, not only these series but also their partial sums have significant meanings in applications (see 
[H[3]). These partial sums are called the multiple harmonic sums (MHS for short) which generalize 
the notion of harmonic sums. In general it is defined as 

JT(s;n):= Yl k i S1 ' ' ' V' > n e Z> . (1) 

l<ki<— <ki<n 

where |s| := Si + • • • + si is called the weight and I = l(s) the length. By convention we set H(s; r) = 
for r = 0, . . . , I — 1, and i?(0; 0) = 1. To save space, we denote by {ei, . . . , e t } d the string formed by 
repeating {e\, . . . , e t ) d times and l d = {l} d - We also write H(s) for the set of all the partial sums 
H (s; n) when no confusion is likely to arise . 

This paper can be considered as a sequel to [21j whose main goal is to provide generalizations of 
Wolstenholme's Theorem for the MHS. Hoffman [H] obtained similar types of results independently. 
One of the results in [21] is the following generalization to homogeneous MHS: 
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Theorem 1.1. [2 1 1 Theorem 2.13]) Let s and I be two positive integers. Letp be an odd prime such 
that p > I + 2 and p — 1 divides none of ks and ks + 1 for k = 1, . . . , I. Then 



H({s} l ;p-1) 



(mod p ) if Is is even, 
(mod p 2 ) if Is is odd. 



In particular, the above is always true if p > Is + 3. 

One can also investigate the sums H(s; n) with fixed s but varying n. Such a study for harmonic 
series was initiated systematically by Eswarathasan and Levine [5] and Boyd [I], independently. It 
turns out that to obtain precise information one has to study Wolstenholme type congruences in 
some detail and so these two directions of research are interwoven into each other rather tightly. To 
state our main results and conjectures we define 

a ( s ' rii 

H(s;n) = ' r , a(s; n), b(s; n) E N, gcd(a(s; n), b(s; n)) = 1. 
o(s; n) 

For completeness, we set a(s; 0) = and b(s; 0) = 1. Fixing a prime p we are interested in the 
p-divisibility of the integers a(s; n) and b(s; n) for varying n. Thus we put i? (s; n) inside Q p , the 
fractional field of the p-adic integers and let v p be the discrete valuation on Q p such that v p (p) = 1. 
In this general situation we're forced to change the notation used by the previous authors. For any 
m £ N and seN 1 , put 

I(s\m) :={n £ Z>o : b(s;n) ^ (mod m)}, 
J(s|m) :={n e Z> : a(s; n) = (mod m)}. 

Note that J(s\m) ^ since £ J(s|m) always. For any prime p we call J(s|p) the p-divisible set of 
the MHS #(s) defined by Eq. ©. 

In [1] Boyd presented a heuristic argument by modeling on simple branching processes to convince 
us that the p-divisible set of the harmonic series is finite for every prime p (this is also independently 
conjectured by Eswarathasan and Levine [H Conjecture A]). Boyd also proves this conjecture for all 
primes less than 550 except for 83, 127 and 397. We now provide a generalization: 

Conjecture 1.2. For any s and any prime p the p-divisible set J(s\p) is finite. 

Although we are not able to prove this conjecture in general, we obtain a lot of partial results. 
The primary tool to prove these when l(s) > 2 is our Criterion Theorem 12.11 Fixing an arbitrary 
prime p we define 

Go = {0} and G t = {n : p^ 1 <n<p*} for t £ N. 

Criterion Theorem. Let I > 2 be a positive integer and p be a prime such that I £ Gt ■ Let 
s = (si, • • • , si) £ N l and put m = min{sj : 1 < i < I}. For t £ N set f(s,p; t) = mm{— v p [H(s; n)) : 
n G Gt}. If there is t > to such that 

f(s,p; t) > (|s| - m)(r - 1) - m, 

then J(s\p) is finite. 

We list only some results obtained by applying our Criterion Theorem below. More examples 
including those when / = 1 can be found in Sec. [2]|4]and in the online supplement [20] . 

Theorem 1.3. The p-divisible set J(s\p) is finite if 

(1) . s = (1, 1) and p = 3,7, 13, 31, or s = (1, 1, 1) and p = 3, or s = (1, 1, 2) and p = 7. 

(2) . s = (4, 3, 5), (5, 3, 4) and p = 17. 

(3) . s = {s} 1 , 1 < I < 20, s > 2, andp = 2. 

(4) . s = (s, t), s, t < 20, t > 2, and p = 2,3, 5. 

(5) . s = (r, s, t), r,s,t < 10, t > 2, and p = 2,3, 5. 
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(6). s = (q, r, s, t), q,r, s,t < 4, t > 2, and p = 2,3, 5. 

Moreover, for s in the last four cases we have J(s|2) = {0}. 

Conjecture 1.4. For every s the 2-divisible set J(s|2) = {0}. 

In [8j Conjecture B] Eswarathasan and Levine state that there should be infinitely many primes 
p (so called harmonic primes) such that J(l\p) = {0,p— l,p 2 —p,p 2 — 1}- Boyd [4 further suggest 
1/e as the expected density of such primes. For any s we extend this notion to define the reserved 
(divisibility) set RJ(s; x) of polynomials in x with rational coefficients. For any prime p > |s| + 3 we 
have RJ(s;p) C J(s\p) and there are primes p (called reserved primes for s) such that equality holds. 
We determine RJ(s) for many types of s in Theorem 17.21 Furthermore, we argue by heuristics that 
the following conjecture should be true. 

Conjecture 1.5. If l(s) = 1 and s — s > 2 then the proportion of primes p with J(s\p) = RJ(s) is 
\j \fe. This proportion is equal to 1/e for all other s. 

This conjecture is supported by very strong numerical and theoretical evidence which we gather in 
[20l Appendix II] and in Theorem 17.21 It also generalizes Boyd's density conjecture of the harmonic 
primes. 

At the end of this paper we put forward some more conjectures of J(s\p) related to the distribution 
of irregular primes. 

This work was partially supported by NSF grant DMS0348258 and the Faculty Development 
Fund from Eckerd College. 

2 A process to determine J(s\p) 

For any positive integer n let n = pn + r, where n, r G N and < r < p — 1. For any seN' define 

l<fci<— <k t <n 1 ' 
p\k 1 k 2 ---k l 

Then by a straightforward computation using the shuffle trick we have: for any s,leN 

H{s-n)=H*{s-n)+p- s • H{s;n), (2) 
/ 

H({s} l ;n) = J2p~ ks ■ H({s} k -h) ■ H*({s} l - k ;n). (3) 

k=0 

where H({s}°; m) = H*({s}°; m) = 1 for any integer m. We omit the proofs of these formulas whose 
main ingredient is contained in the proof of the main Criterion Theorem 12.11 below. Both of these 
formulas are generalizations of [51 (2.2)] for partial sums of harmonic series. They are the primary 
tools to study Conjecture 11.21 for homogeneous MHS. 

For more general MHS we need a more complicated version of these formulas. Fixing an arbitrary 
prime p we define 

G = {0} and G t = {n : p*" 1 < n < p 1 } for t G N. 
For any s e N 1 we set J t (s|p) = G t H J(s\p). 

Theorem 2.1. (Criterion Theorem) Let I > 2 be a positive integer and p be a prime such that 
I G G to . Let s = (si,-- - ,Si) € N' and put m = min{si : 1 < i < I}. For t £ N set f(s,p;t) = 
min{— v p (H (s; n)j : n G Gt}. If there is t > to such that 

/(s,p;r) > (|s| -m)(r- 1) - to, 

then J(s\p) is finite. 
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Proof. Let ri = pn |r£ G T +±. By definition we have 



H(s;n)= EE E gf x.!. gr - ( 4 ) 



a+/3=l l<i 1 <--<i a <n l<Ki<-<Ki<n 

a,/3>0 l<jl<-<3i3<n p\K ia ,l<a<a 

{i 1 ,...,i a }U{j 1 ,...,j }={l,...,n} p\K 3b ,l<b<f3 



Note that the terms corresponding to a = I, (3 — form the series A — H(s; n) /p' s ^ since r > to 
so that n > p*° > i. For all other terms with j3 > 1 we have the natural bound 

« p (^"--^ S! )<(N-m)r 

since Kj < p T+1 for all j and one of them is prime to p. Set B = H(s; n) — A then obviously 
v p{B) > — (|s| — m)r. Since n G G r by assumption of t we know that 

v p (A) < -|s| - (|s| - m)(r - 1) +m = -(|s| - m)r < u p (B). 

By induction it's easy to see that for alH > r and n € Gt we have 

v p (H(s;n)) < -\s\(t - r) - (|s| - m)(r - 1) + m = -|s|(< - 1) + tot < 0. 

This shows clearly that J(s\p) is finite. □ 

Corollary 2.2. Let s and Z be two positive integers, I > 2, and p be a prime. Suppose I € Gf . T/ien 
the p- divisible set J({s} l \p) is finite if there existsr > to such that fi(r) := f{{s} l ,p;T) > (1—1)st—s. 

Proposition 2.3. Let r, s, t be three positive integers. Let s = (r, s) with 1 < r < 10 and 2 < s < 10. 

ors = (r,s,t) with 2 < r, s, t < 5. Then there's always some prime p > \s\ +3 such that the p- divisible 
set J(s\p) = RJ(s;p) is finite where RJ(s;p) is given in Theorem 17.21 

Proof. The set RJ(s;p) will be defined for general s in Definition [7j] and computed in Theorem l7.2t 
The proof of the proposition follows from the Criterion Theorem 12. II by computation. To save space 
we put the details online [2U]. □ 

3 Finiteness of J(s\p) 

We now describe an approach to determine the p-divisible set J(s\p) for any given positive integer 
s and odd prime p. This is essentially discovered by Eswarathasan and Levine [5] and by Boyd [3] , 
independently. It follows quickly from Eq. that 

n € I(s\p) if and only if n £ J(s\p s ). (5) 

n £ J(s\p) implies h £ J(s\p s ). (6) 
I(s\p) = pJ(s\ P s ) +R, R = {0, 1, • • • ,p - 1}. 



Therefore 

It's also clear that 



Remark 3.1. The case when s > 1 is very different from that of s = 1 considered by previous 
authors in that the information of I(s\p) is in general not enough to determine J(s\p). 

To get an equivalent condition of Eq. ^ we need a partial generalization of [8, Lemma 3.1]. Set 
the parity function p(m) = 1 if to is odd and p(m) = 2 if m is even. 

Lemma 3.2. Let p be an odd prime and s a positive integers. If p — 1 \ s, s + 1 then we have 

H*{s;pn) = (mod p^-^). (7) 
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Proof. By definition 

„ 1 2zJ / „ 1 \ 

(8) 



l<fe<n m=0 \ mp<k<(m+l)p 



(p,k) = l 

The lemma follows from the fact that each inner sum in the parentheses satisfies the congruence in 
Eq. J7]) which can be proved by the same argument as that in the proof of |211 Lemma 2.2] (when n 
is odd the shorter proof suffices). It also follows from Corollary 1]. □ 

Proposition 3.3. Let s,IeN and p be an odd prime such that p > I + 2 and p — 1 divides none of 
ks and ks + 1 /or fc = 1, .... I. Then 

H*({s} l ;pn) = (mod p^" 1 )). 

The proof as well as the result itself is similar to that of Theorem 1 1.1 1 so we leave the details to the 
interested reader. The first step of induction is given as Lemma 13.21 above. In fact, the proposition 
itself reduces to Theorem 11.11 when n = 1 . 

Definition 3.4. For n £ J(s\p s ) there is a unique integer ip s (s\p;n) 6 [0,p— 1] such that 

ip s (s\p;n) = — H(s;n) (mod p). (9) 

Lemma 3.5. For n = pn + r, < r < p, we have 

H*(s;n}-H*(s;pn)=H(s;r) (mod p). (10) 

Furthermore, if n € J(s|p s ) i/ien 

if (s; n) = i? (s; r) + tp s (s\p; n) (mod p) (11) 

Proof. This follows from Eq. ((2|) and Lemma [3~2l Also see the proof of [8j Lemma 3.2]. □ 

Theorem 3.6. Let n = pn + r, < r < p. Then n € J(s\p) if and only if 

n e J(s\p s ) and H(s;r) + tp s (n) = (mod p). (12) 

Proof. If n E J(s\p) then Eq. ([6|) implies that n S J(s|p s ). In addition, congruence (fT2|) follows 
immediately from Eq. (fTT|) . On the other hand, if Eq. (fT2|) holds then Eq. (fTT|) implies that n € J(s\p) 
and the proof is complete. □ 

We now use the above theorem to define a branching process by using the sets Gt which will 
compute J{s\p) if it's finite. 

Proposition 3.7. Let s be a positive integer and p an odd prime. Then J(s\p) = U^ Jt(s|p) where 
Jt(s\p) can be determined recursively by 

Jt+i(s\p) = {n = ph + r :h€ J t (s\p s ),r G R,v p (H(s;r) + ip s (n)) > 0} 

for f 6 N. Here, as before, R = {0, 1, • • • ,p — 1}. 

The next corollary follows naturally. 

Corollary 3.8. Let s be a positive integer and p an odd prime. Then J(s\p) is finite if and only if 
J t (s\p s ) = for some t € N. 



An easy computation according to Corollary 13.81 yields the following concrete result. 
Proposition 3.9. Let s be a positive integer. Then J(s\p) is finite for primes p = 2, 3, 5, 7. 
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Proof. (1) p = 2. We claim that J(s\2) — {0}. We can prove that 2 does not divide H(s;n) by 
induction on n. This is clear for n = 1 and n = 2 because H(s;l) — l,H(s;2) = (1 + 2 S )/2 S . 
Suppose r ^ ^(s|2) for all r < n and n £ J(s\2). If n is odd then let H (s; n — 1) = ^ where a is odd 
by inductive assumption. Then 

„. . a 1 Na + 2B 
H (s; n) 



2b n s l.c.m(26,n s ) ' 

where N — n s /gcd(n s , b) and B = &/gcd(n s , 6). Hence Na + 2B is odd because both N and a are 
odd, which is a contradiction. If n — 2n then 

" / 1 1 \ 1 

2 -) - £ ( jtk^w + PF J s fi + ^ i/(s; " } (mod 2) ' 

By inductive assumption 2 { _ff (s; n) which implies that 2 { (s; 2n). So n can not belong to J(s\2) 
either if n is even. This shows that J(s\2) = {0}. In fact, it is not hard to see that for n £ Gt, t>l 
we have 

V2(H(a;n))=-(t-l)a. (13) 

For 3 < p < 7 Eswarathasan and Levine [S] have shown that J(l\p) are finite. We also know that 
when s < 4 then J(s\p) are finite for these primes by explicit computation [20]. Assume s > 4. Then 
by Corollarv l3.8l we only need to show that Ji(s\p s ) = 0. We need [2 1 1 Corollary 2.7] which implies 
that if p > 3 is a regular prime then 

H(s;p- 1)^0 (mod p s ) for s > 4. (14) 

(2) p = 3. Neither H(s; 1) = 1 nor H(s; 2) = 1 + 1/2 S is divisible by 3 s . so Ji(s|3 s ) = 0. 

(3) p = 5. Neither iJ(s; 1) = 1 nor F(s; 2) = 1 + 1/2 S is divisible by 5 s . Now 



6 s H(s; 3) = 2 s + 3 s + 6 s = 2 s + (-2) f 



2 • 4" + 1 (mod 5) if s = 2n, 
1 (mod 5) if s is odd. 



So we always have H(s; 3) ^ (mod 5), i.e., 3 g Ji(s|5 s ). Finally, Eq. ((H} implies that 4 g Ji(s|5 s ) 
for s > 4 because 5 is a regular prime. Hence Ji(s\5 s ) = 0. 

(4) p=7. Clearly 1,2 ^ Ji(s\7 s ) and 6 s i7(s; 3) = 2 s + 3 s + 6 s < 7 s when s > 4. Now 

ff( S ; 4) = H(s; 6) - 1 - 1 = (-l) s+1 (l + ^ (mod 7). 



Because 2 3 = 1 (mod 7) we get 

1 



(mod 7) if s = (mod 3), 
1 + — = <^ 3/2 (mod 7) if s = 1 (mod 3), 
(mod 7) if s = — 1 (mod 3). 



Therefore 4 g Ji(s|7 s ). Similarly, ff(s; 5) = (-l) s+1 ^ (mod 7). Finally, it follows from Eq. dHJ) 
that 6 g Ji(s|7 s ) for s > 4. These show that Ji(s|7 s ) = for all s > 4. □ 

Remark 3.10. The case p = 11 is not so easy since if (3; 4) = (mod 11) and moreover, for any 
positive integer e there is some s < p e (p — 1) such that H(s; 4) = (mod ll e+1 ). 

We also computed J(s\p) for some other s and p (see [ID])) which confirms the following 

Proposition 3.11. Let p be a prime such that p < 3001. Then J(s\p) is finite for 2 < s < 300. 
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4 Finiteness of J({s} l \p) 

In order to apply Criterion Theorem 12.11 we set 

fl(t) := f({s} l ,p;t) = min{-v p (H({s} l ;n)) : n e G t } Vt > 1. 
We first look at the case s > 2. 

Lemma 4.1. For all s > 2 we have «2(3 S + 1) = p(s — 1) which is 1 if s is even and 2 if s is odd. 
In particular, we always have 3 s + 1 ^ (mod 2 s ). 



Proof. This is clear because 

3 s + 1 = 



9™ + 1 = 2 (mod 8) if s = 2n, 

3-9" + 1 = 4 (mod 8) ifs = 2n+l. 



□ 



Proposition 4.2. Let s > 2 and I < 20 be two positive integers. Then the 2-divisible set J({s} l \2) 
is finite. 



Proof. When / = 1 this is included in Proposition 



#(s;2) = l + i ff( S ,s;2) = i H(s;3) 



So we assume s, I > 2. Then 
6 s + 3 s + 2 s 



6 s 



Further, by Lemma l4.ll we know that 



1 \ 1 3 s + 1 + 2 s 



G* 



has at least a factor 2 in the denominator. Therefore we can take t — 2 to get /(r) > 1 > 
(2s — s)(r — 1) — s = 0. So the condition in Corollary |2.2l is satisfied and consequently J(s, s\2) = {0}. 

A detailed study of using Lemma |4~T1 tells more. Let t > and n G Gt+2- Then by induction and 
Eq. we can easily show that 



v 2 (H(s, s; n)) 
Putting I = 3 in the following equation 



-(2t+l)s if 2-2* < n < 3- 2*, 

p(s- 1) - (2t+ l)s if 3-2* <n<4-2*. 



H({ S } 1 ; n) = H({s} l ;n - 1) + — ^({s}'- 1 ; n - 1), 



and applying induction on t we can show that 
v 2 (H(s,s,s;n)) 



p(s-l)-3ts if4<2n/2*<5, 
-3ts if 5 < 2n/2* < 6, 

-(3i + l)s if6<2n/2*<8. 



(15) 



(16) 



(17) 



So we get J(s, s, s\2) = {0} when s > 2. 

When / > 4 we can utilize Eq. (Tf6|) again. However, even in the case I = 4 it is very complicated 
already. Nevertheless the idea is straightforward so we omit the details of the proof. Suppose s = 2 
and n £ Gt+2 with t > 1 (note that H({s} 4 ; n) = for all n < 3). Then we have 



U2 (F({2}V)) = 



f(l) 


- 2(4t 


-1) 


if 32 < 16n/2* < 48, 


(2) 


- 8t 




if 48 < 1671/2* < 56, 


(3) 


-8t + 


S(t) 


if 56 < 1671/2* < 57, 


(4) 


-8t + 


7 


if 57 < 16n/2* < 58, 


(5) 


- 2(4* 


-2) 


if 58 < 16n/2* < 60, 


,(6) 


- 2(4* 


-1) 


if 60 < 16n/2* < 64. 



(18) 
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Here if * = 1 then (3)- (6) merge into (6); if * = 2 then S(t) = 5 and (3)- (5) merge into (3); if t = 3 
then (3) and (4) merge into (3); if t > 3 then 5(t) = 6. When s = 3 and ra e Gt+2 with t > 1 we have 



(1) 


- 3(4* - 


1) 


if 8 < 4n/2* < 12, 


(2) 


- 12* 




if 12 < 4n/2* < 14, 


(3) 


- 3(4* - 


1) 


if 14 < 4n/2* < 15, 


(4) 


- 3(4* - 


+ 1 


if 15 < 4n/2* < 16. 



Here if * = 1 then (3) and (4) merge into (4). When s > 4 we have 
v 2 (H({ S y ; n)) - 



(19) 



s(4* - 1) 
4s* 
4st + 2p(s 



if 4 < 2n/2* < 6, 
if 6 < 2n/2* < 7, 
if 7 < 2n/2* < 8. 



(20) 



Equations (fl8 ]) -([2D |l imply that J({s} 4 \2) = {0} for all s > 2. 

Similar computation shows that when I — 5 and n € Gt+2 with * > 1 we have 



v 2 (H({2} 5 ;n)) 



(1) 
(2) 
(3) 
(4) 
(5) 
1(6) 



- s(5* - 3) 

- s(5*-3) 

- s(5*-3) 
-s(5*-2) 
-s(5*-l) 
-s(5*-l) 



2p(s-l) 
3 



1 



if 16 < 8n/2* < 17, 
if 17 < 8n/2* < 18, 
if 18 < 8n/2* < 20, 
if 20 < 871/2* < 24, 
if 24 < 8n/2* < 28, 
if 28 < 8n/2 t < 32. 



(21) 



Here if t = 1 then (l)-(3) do not appear; if s > 3 then (1) and (2)) merge into (1). This implies that 
J({s} 5 |2) = {0} for all s > 2. 

As the length becomes longer (i.e. I gets bigger) there are more and more cases. The number 
of cases, denoted by C(l), is independent of s when s is large enough and tends to increase with / 
though not always. We compute the following 

C(6) = 5, C(7) = 7, C(8) - 6, C(9) = 8, C(10) = 8, 
C(ll) = 11, C(12) = 10, C(13) = 12, C(15) = 15, C(16) = 12, 
C(17) = 15, C(18) = 14, C(19) = 18, C(20) = 15. 



After tedious verification we find J({s}'|2) = {0} for all I < 20 and s > 2. 



□ 



For any given / by similar method we should be able to determine J({s} 1 12) for all s > 2. However, 
for odd primes p we can only extend this result to J({s} l \p) for small I and small s with the aid of 
computers. 

Proposition 4.3. Let s and I be two positive integers. Suppose 2 < s < 10 and 2 < I < 10. Then the 
p-divisible set J({s} l \p) is finite for the consecutive five primes immediately after Is + 2. Moreover 
there's always some prime p such that J({s} l \p) = RJ({s} l ;p) where 

RJ({s} l :p) = 



{0,p-l} if2\s, 
{0,i+(p-l)/2,p-l:0<i<l-l} if2\s. 



Proof. The set RJ(s;p) will be defined for general s in Definition [7j] and computed in Theorem [721 
The proof of the proposition follows from Corollary 12. 21 by computer computation. To save space we 
put the details online [5D]. □ 

In the rest of this section we turn to the case s = 1. We may assume I > 2 since the harmonic 
series has been handled by [H] and [J] . According Corollary 12.21 if we can find r large enough such 
that fi(r) >(l- 1)(t - 1) then J(l l \p) is finite. 
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Proposition 4.4. 1. The p-divisible set J(l 2 \p) is finite if p — 3, 7, 13, 31. 

2. Let s,t < 20 and t > 2. TTien i/ie se£ J(s, t|f>) is finite for p — 2, 3, 5. 

5. Lef r, s,t < 10 and t > 2. Then the set J(r, s, t\p) is finite for p — 2, 3, 5. 

^. Let q,r, s,t < 4 and t > 2. Then the set J(g, r, s, t|p) is finite for p = 2, 3, 5. 

Proof. We only need to find r satisfying the condition of Corollary 12. 21 
(1) For each r in the following we have /2(t) = t — 1. 

p = 3. Take r = 6. Then computation reveals that J(l, 1|3) = {0, 5}. If I = 3 then we take r = 10. 
Then we have / 3 (r) > 2(r - 1). Note that in G 9 there is n = 17770 such that v 3 (H(l 3 ;n)) = -15 
so /a(9) = 15. By Corollary 12. 21 and simple computation we see that J(l 3 |3) = {0,8}. 

p=7. Take r = 4. Then J(l,l|7) = {0,4,6,7,13}. 

p = 13. Take r = 4. Then J(l, 1|13) = {0, 12, 13, 25}. 

p = 31. Take r = 4. Then J(l,l|31) = {0,17,22,30,31,61}. 

For the last three cases with p — 2, 3, 5 we put the result of computation online [5D] . For example, 
we can take r = 10 and show that J(l, 1, 1|3) = {0,8}. □ 

Remark 4.5. We could extend our results to larger I and some other primes p but it would be very 
time consuming with our slow PCs. However, even in the case s = (1, 1) similar process fails for 
p = 2. Computations suggest that J(l, 1|2) = {0}, J(l,l|5) = {0,4,5,9}, J(1,1|H) = {0,10,11,21} 
and J(l, 1 1 IT) = {0, 11, 13, 16, 17,33}. We will analyze the situation for p — 2 in detail in the next 
section. 



5 Sequences related to J(s, 1|2) 



One may wonder what goes wrong in Proposition [4T] if we let s = (1,1) and p — 2. We will see that, 
amazingly, this problem is closely related to some pseudo-random process. 

Only in this section we adopt the shorthand Hi(n) := H(l;n) and H2{n) := H(l,l;n). Let's 
start with the first few partial sums of H2(n) when 2 < n < 14. Here ~ means we only consider the 
fractional part of the numbers. 

1 11 7 23 

ff 2 (2)~-, # 2 (3) ~1, fTa(4)~— , ff 2 (5)~-, H 2 (6) ~ 

H 2 (7)^^, #2(8)-^-, F 2 (9)~i^, #2(10)~™ 
2K ' 180' 1 ; 10080' v 7 2016' 1 ; 50400' 

ff2 (ll)~^, #2(12) ~ ff a (13)~i^, # 2 (14)~^. 

v 7 50400 v ; 50400 v ' 415800 v ; 14850 

It looks like 2 never divides the numerator and moreover, the 2-powers in the denominators of Hi(n) 
tend to increase with n, though not always. To proceed we need to know the 2-divisibility of Hl{n). 



Lemma 5.1. Let n be a positive integer. Then 

H* x (n) 



(mod 4) if n = 0, 3 (mod 4), 

1 (mod 4) i/n = l,2 (mod 4). 



Proof. If n is even then obviously H*(n) — H*(n— 1). So we only need to consider n = 1, 3 (mod 4). 
Set 8 = 1 if n = 41 + 1 and 5 = if n = 41 - 1. Then 

H *r n) = J?_ + y_!_ = J?_ + Y(_}_ + l )= 6 (mod 4) 

iy ' 41 + 1 ^ 2i - 1 4Z + l ^V2i-14Z-2i + l/ v y 

i=l i—l v 7 

as desired. □ 
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Remark 5.2. By working more carefully we can obtain the following improvement of Lemma [5.11 
if n = 2 d m or n = 2 d m — 1 where m is odd and d > 1. Then v 2 (H*(n)) = 2(d — 1). However, the 
proof is complicated and it is not needed in the rest of the paper so we leave the proof of this general 
statement to the interested reader. 

The following result is exactly the reason why Corollary 12. 21 cannot be applied to J(l, 1|2). 

Proposition 5.3. For any t > 2, there is a unique n t <E G t such that v 2 (H 2 (n t )) > 2 — i whereas 
for all n t ^ n € Gt we have v 2 (H 2 (n)) < 1 — t. Therefore, for all positive integers n ^ {n t } t >i the 
numerator of H 2 (n) is an odd integer. 

Proof. Note that G\ — {1} and G2 — {2,3}. Thus n 2 = 3 because -ff 2 (3) = 1. Assume that t > 3 
and each ni has been found in Gi uniquely for i < t. Let n — 2h + r G Gt+i for r = or 1. When 
I = p = 2 and s = 1 Eq. J3]) becomes 

H 2 (n) = H* 2 (n) + ^H 1 (h)H^(n) + ~H 2 (n). (22) 

It's easy to show that V2{Hi(m)) = 1 — t for m € Gt by induction and the recursive relation 
Hi(n) = H*(n) + Hi(h)/2. If n ^ n t then we have V2{H2{nj) < 1 — i and hence 

«2(fl2(n)) = min{«a(ff?(n)) - 1, -1 - 1} = -1 - 1 < 1 - (t + 1). 

Suppose now h = n t and n = 2n t + r t . We consider four possible cases. 

(i) If V2(H 2 (n t )) = 2 — t and ^(-/^(n)) > 1 then n ^ n t+1 because 

v 2 {H 2 {n)) = mm{v 2 (H*(n)) - t, -t} = -t = 1 - (t + 1). 

(ii) If v 2 {H 2 {nt)) = 2 - t and ^(^(n)) = then 



Hence n = n f+ i because 



« 2 (ff 2 (n)) > 1 - t = 2 - (* + 1). 



(hi) If V2 {H2{ntj) > 3—ta.ndv2(Hi(n)) = Othenn 7^ because v 2 (H 2 (n)) = min{u2(.ff*(n))- 
i, «2 (fT 2 (nt)) - 2 - 1} = -t = 1 - (t + 1). 

(iv) If v 2 (H 2 (nt)) > 3 — t and U2 (if* (n)) > 1 then n = n f+ i because 

v 2 (H 2 {n)) > mm{v 2 (H*{n)) - t,v 2 (H 2 (n t )) - 2 - t} > 1 - t = 2 - (t + 1). 



Now if n t = 2Z is even then by Lemma 15.11 

(1) 2n t + 1 = 1 (mod 4) and w 2 (iJ*(2n t + 1)) = 0, and 

(2) 2n t = (mod 4) and v 2 (H*(2n t + 1)) > 1. 
If n t — 21 + 1 is odd then by Lemma [5~T1 

(3) 2n t + 1 = 3 (mod 4) and w 2 (-ff*(2n t + 1)) > 1 and 

(4) 2n 4 = 2 (mod 4) and v 2 (H$(2n t + 1)) = 0. 
Therefore, we have four situations to consider: 

(a) n t is even and v 2 (H 2 (nt)) = 2 — t. Then n t +i = 2n t + 1 by (1) and (ii). 

(b) nt is even and v 2 (H 2 (nt)) > 3 — t. Then n t +\ = 2n t by (2) and (iv). 

(c) nt is odd and v 2 [H 2 {rit)) > 3 — t. Then n t +\ = 2n t + 1 by (3) and (iv). 

(d) nt is odd and v 2 (H 2 (rit)j = 2 — t. Then n t +\ — 2n t by (4) and (ii). 

It follows that n t +i £ Gt+i is uniquely determined. This finishes the proof of the proposition by 
induction. □ 
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Denote the dyadic valuation u 2 {Hi (nt)) by —Wt- Then we have the following two interesting 
sequences: 

{n t } t >2 = {3, 6, 13, 27, 54, 109, 219, 439, 879, 1759, 3518, 7037, 14075, 28151, 

56303, 112606, 225212, 450424, 900848, 1801696, 3603393, . . . } (23) 
{w t }t>2 = {0, 1, 3, 4, 3, 3, 5, 7, 9, 10, 9, 10, 12, 14, 13, 13, 15, 17, 19, 19, ... } (24) 



Set n = r 2 = 1 and define r t = or 1 for t > 3 as determined in the proof of Proposition 15.31 such 
that rit+i — 2n t + r t . Then clearly nt can be written as 

n t = (rir 2 ■ ■ .r t ) 2 (25) 

in binary system and apparently the sequence {nt} increases very fast. Further, the occurrence of 
r t = or r t — 1 does not seem to have any predictable pattern and in fact we believe it is related 
to some pseudo-random process. By this we mean the following. First we of course conjecture that 
{wt}t>3 is always bounded below by 1 which is equivalent to say J(l,l|2) = {0}. We also have 
proved that Wt is bounded above by t — 1 and it is not hard to see that w± = t — 1 for infinitely many 
t's. It's also conceivable that wt are near t — 1 most of the time. However, we believe Wt could move 
very far away from t—1 for very large t. At the present stage, we could not even determine whether 
the difference between wt and t—1 can be arbitrarily large. 

Remark 5.4. We put the two sequences {nt]t>2 and {wt}t>2 in Sloane's online database of integer 
sequences as A079403(rt) and A079404(n), respectively. Shortly after Benoit Cloitre emailed me a 
formula for the known terms of {n t }t>2- 

n t = L2*- 1 cJ (26) 
where c = 1.718232.... Indeed, it's easy to see that 

ni = l,n t = 2^n(l + g : ). V*>2. 

Further, 



t— >oo J- J- y 2nk I 

k—l v / 

exists by comparison test. From Eq. (|25p we get 

c = (ri.r 2 r 3 r 4 . . . ) 2 = (1. 101101111101111000001. ..) 2 = 1.718232... 
Moreover, using binary system we see that the integral part of 2 t_1 c is exactly n t , as desired. 
We can easily generalize Proposition to the following. 



Proposition 5.5. For any t > 2, there is a unique n t G Gt such that u 2 (iJ(s, 1; n t )) > — s(t— 1) + 1 
whereas for all n t ^ n E G t we have u 2 (i? (s, 1; n)) < —s(t — 1). Therefore, for every positive integer 
n (ji {nt}t>i the numerator of H(s, is not divisible by 2. 

Proof. We can assume that s > 2 because of Proposition 15.31 The key to the proof is Eq. Q which 
yields that 

H(s, 1; n) = H*(s, 1; n) + U(s, 1; n) + V(s, 1; n) + ^H(a, 1; n) 

where 

l<2k<l<n,2\l l<k<pl<n,2\k 

Now it's easy to see that if v 2 (7J(s, 1; nj) < —s(t — 1) then u 2 (iT(s, 1; n)j < — st. 
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When t — 2 we find n 2 — 2 always because 

, 1 , , 1 2 s + 1 
H(s,l;2) = -, H(s,l;3) = - + Y -^ r . 

Assume that t > 3 and n is the unique n t G Gt such that v 2 {H(s, 1; n)) > — s(< — 1). Then 
V2(H(s; ri)) = — s(i — 1) by Eq. (fT3|) . So we can always uniquely choose r t so that for n = 2h + r t € 
Gt+i 

v 2 (U(s,l;n)) = v 2 (U(s,l;2n) + rH(s;n) /2 s ) = -st 

if v 2 (H(s,l;n)) = -s(t -1) + 1. If v 2 (H(s, 1; n)) > -s(i-l) + l then we can uniquely choose r t so 
that 

v 2 (U(s,l;n)) > -st+1. 

The upshot is for nt € Gt there is a unique n t +\ G Gt+i satisfying the condition of the proposition. 
This finishes the proof. □ 

In general, we cannot apply Criterion Theorem to determine the finiteness of J((s, 1)|2) for any 
s £ N 1 , because of the existence of similar sequences. Moreover we believe 2 never divides the 
numerator of any multiple harmonic sum. 

Conjecture 5.6. Let I be an arbitrary positive integer and se N 1 . Then the 2-divisible set J(s|2) = 
{0}. 

We have verified this conjecture for all s = (s, t) and s ~ (r, s, t) with 1 < r, s < 10 and 2 < t < 10, 
and for all s = (u,r,s,t) with 1 < u, r, s < 4 and 2 < t < 4. See [20]. The computation is very 
time-consuming, for example when s = (1,4,4,2) the Maple program runs more than 3.5 hours on 
my PC with Pentium 4 CPU 3.06GHZ and 512 MB RAM. The same program in GP Pari runs a 
little faster. We put the program at the end of our online supplement [20] . 

We believe that among all possible s and prime p the cases (s, 1) are the only ones that our 
Criterion Theorem fails (see [20]). Let me sketch a heuristic argument for this belief for the case 
s = l l and p > 3. 

By the recursive relation 

/ 

H(l l ;n) = ^P _J -H(l j ;n) ■ H* ;n). 
j=o 

it is not hard to see that the size of r we are looking for in the Criterion Theorem depends on the 
length of the sequences {nt}t>t not satisfying the condition in the theorem, where n t € Gt and 
nt+i = pn t + r t for some < r t < p. If n t is already found then the existence of nt+i depends 
on H*(l;n) essentially, which we assume to distribute among (p + l)/2 values modulo p by the 
symmetric structure of J\{l\p) (see Sec. O. So nt produces two possible nt+i or no n t +\ with the 
same probability q = (p — l)/(2p), and it produces exactly one n t +\ with probability 1/p. 

Let's assume that a certain cell reproduces itself according a similar law as above, namely, it 
clones itself or dies in the next generation with the same probability q, and with probability 1/p it 
stays alive without reproduction. Let pk be the probability that starting from fc cells in the beginning 
the cells eventually all die out. We claim that pk = 1 for all k. Indeed, it is not too hard to see that 
we only need to showpi = 1. This follows from the Criticality Theorem for Galton- Watson branching 
process (see [TUl Preface] or pQ p. 7, Theorem 1]) because the average offspring is 2q+ 1/p = 1. 

6 The structure of Ji({s} l \p) 

Set Jt(s\m) = {0} U J t (s|m) for any positive integers t, m and s e N z . The next result is easy 
but very useful in determining the structure of Ji(s\p) since it tells us essentially that J°(s|p) is 
symmetric about (p — l)/2. 

Proposition 6.1. Let p be an odd prime and s € N . Let r € {1, . . . ,p — 2}. Then r € J\{s\p) if and 
only if p — 1 — r € Ji(s\p). 
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Proof, lip — l\s then Ji{s\p) — because H(s; r) = r (mod p) for all r £ {1, . . . ,p — 2}. If p — 1 f s 
then we have 

k=l k=p-r v ^ ' k=p-r 

Subtracting = (— l) s X)fe=i w ( m °d P) from the above we get the desired result. □ 

Remark 6.2. We feel prompted to mention that the symmetry of J°(s|p) is not enjoyed by Jj 3 (s|p 2 ). 
For instance, while J x (5|37) = {0, 6, 9, 12, 18, 24, 27, 30, 36} is symmetric about 18 the set J!°(5|37 2 ) = 
{0, 6, 36} is not. 

Now that we know Jf(s\p) is symmetric we may wonder what happens to the center (p — l)/2. 
When s is odd the answer is related to the irregularity of primes. 

Proposition 6.3. Let p be an odd prime and s be a positive integer such that p — 1 \ s. Let n be the 

unique positive integer such that s = n (mod p — 1) and 2 < n < p — 2. // n < p — 3 then we have 

!^=^- (modp), ifn = l; 

^B p _ n (mod p), ifn>lisodd; (27) 

" ( 2(«+i) ] P B p-n-i (mod p 2 ), if n is even. 

Therefore, 

(a) If s is even then (p — l)/2 £ Ji(s\p). 

(b) If s is odd and (p,p — n) is an irregular pair then (p — l)/2 £ Ji(s\p). 

(c) If s is odd and (p — l)/2 £ Ji(s\p) and 2 s ^ 2 (mod p) then (p,p — n) is an irregular pair. 
In particular, if s > 3 is odd and p > 2 s — 2 then (p — l)/2 £ Ji(s\p) if and only if (p,p — s) is an 
irregular pair. 

Proof. The congruence (|27|) is essentially [17l Cor. 5.2(b)]. The rest follows immediately. □ 

Remark 6.4. Proposition 16.31 tells us that when s = 1 (mod p) then (p — l)/2 £ J(s\p) if and only 
if p is a Wieferich prime. There are only two known such primes: p — 1093 and p = 3511. If any 
others exist, they must be greater than 1.2510 15 according to [15] . 

The above proposition says that if s is even and p — 1 \ s then (p — l)/2 £ J\{s\p). A natural 
question is when (p — l)/2 £ Ji(s\p 2 )l The answer is given below. 

Corollary 6.5. Suppose s is a positive integer such that p ~ 1 j s. Suppose n is the integer between 
2 and p — 1 such that n = s (mod p — 1) and n < p — 4. For n even H(s; (p — 1) /2) = (mod p 2 ) if 
and only if either (p,p — n — 1) is an irregular pair or p\2 n+l — 1. If n odd, p\2 n — 2 and (p,p — n) 
is an irregular pair then H(s; {p — l)/2) = (mod p 2 ). 

Proof. This is clear for n even. If n is odd then the corollary follows from [T71 Thm. 5.2(b)]. □ 

Remark 6.6. For every positive even integer s and every irregular prime p > s + 4 up to 100,000, 
p 2 always exactly divides H(s; (p— l)/2). Is this true in general? The answer is no. A calculation by 
Maple shows that for the 5952nd irregular pair (p,p— n— 1) = (130811, 52324) we have n = 78486 and 
2" +1 = 2 (mod p) and therefore p 2 \H(n+l; (p-l)/2) &ndp 3 \H{n; (p-l)/2) by [10 Thm. 5.2(a), (b)]. 
The peculiarity of this pair was already noticed in [6]. The next two such pairs are (599479,359568) 
(see [7]), and (2010401,1234960) (see [5]). Note that apparently this problem is not related to the 
problem of 2 P = 2 (mod p 2 ). 

Theorem 6.7. Let s be a positive integer and p > 2ls +1 be an odd prime. Then 
{p-l,j+(p-l)/2:j = 0,l,---,l-l}c Ji({2s} l \p). 
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Proof. Let m — (p — l)/2. By [211 Lemma 2.12] there are integers c\ such that 

UH({s} l ;n)= c xH x (s;n), (28) 

\eP(i) 

where P(l) is the set of partitions of I, -ff(Ai,....A r )( s ; n ) — IIj=i H(Xjs; n) and C(;) = (— l) z 1 (Z — 1)!. 
Plugging in n = m we get m £ Ji({2s}'|p) by Proposition 16.31 By definition iJ({s}';g) = for 
q = 1, • • • ,1-1. When g = 1 this implies that ^Aepm c a = by Eq. 1(55)). Hence m+1 £ Ji({2s}'|p). 
Similarly, because 1, l/2 2s , • • • , 1/(1 — l) 2s are linearly independent when regarded as a function of 
s, we see that for all independent variables x±, . . . , Xj, j < I — 1, we have 

E '-aII^ • •••• r . x : 

X=(A 1 ,...,A r )eP(d) 3=1 

The theorem now follows from setting Xj = l/(m + j) 2s for j = 1, • • • , I — 1. □ 

Corollary 6.8. Let s be a positive integer and p > 4s + 1 be an odd prime. Then 

(1) If s is even then (p - l)/2, (p + l)/2 S Ji(s,s|p). 

(2) If s is odd and (p,p— s) is an irregular pair then (p — l)/2 € Ji(s,s|p). Further, if s is odd, 
2 s ^ 2 (modp), anc? (p— l)/2 G Ji(s,s|p), t/ien (p,p — s) is an irregular pair. In particular, if s > 3 
is odd cm<i p > 2 s — 2 t/ien (p — l)/2 € Ji(s, s|p) i/ and only if (p,p — s) is an irregular pair. 

Proof. Let s and p be the integers satisfying the conditions of the corollary. When s is even the 
corollary follows from Theorem 16.71 If s is odd then by Proposition 16.31 and the shuffle relation we 
have 

2H(s, s; (p - l)/2) = H(s; (p - l)/2) 2 - H(2s; (p - l)/2). (29) 
The corollary follows immediately. □ 

7 Reserved set of MHS 

In Conjecture B of [8] Eswarathasan and Levine state that there should be infinitely many primes 
p such that the divisible set J(l\p) = {0,p — l,p 2 — p,p 2 — 1}. Boyd [4j further suggest 1/e as the 
expected density of such primes. The most important steps are to elucidate the structure of Ji(l|p) 
and determine the relation between Jt(l|p) and J t+ i(l|p) for t > 0. We put forward some similar 
results and conjectures concerning the divisible sets of general MHS in this last section. 

Definition 7.1. For any s 6 N z there are finitely many function fo(x) = 0, fi(x), . . . , f r (x) G Q[x] 
such that for all primes p > wt(s) + 3 

/o(p) < fx(p) < ■ ■ ■ < fr(p) and f a (p), fi(p), f r (p) S J(s\p). 

We call the largest r the reserved (divisibility) number of MHS H(s), denoted by p(s). We call the 
corresponding set {fo(x), . . . , / p ( s ) (x)} the reserved (divisibility) set of C( s )j denoted by RJ(s) = 
RJ(s; x). Its t-th segment is RJ t (s) = {f(x) € R J(s) : f(p) < p* — 1 for all prime p} for t>l. Note 
that G RJ t (s) for all t > 0. If J(s\p) = RJ(s;p) for some prime p then is called a reserved prime 
for MHS H(s). 

For example, the reserved number of the harmonic series is 3, the reserved set is R J(1) = {0, x — 
l,x 2 — x,x 2 — 1}, and 5 is a reserved prime for the harmonic series because J(l|5) = {0,4,20,24}. 
We summarize all known reserved sets in the following theorem. 

Theorem 7.2. Let r,s,t,l < 5 be positive integers. Then 

(1) . RJ(1) = {0,x- l,x 2 - l,x 2 -x}. 

(2) . If 1 + 1 is odd then RJ w (l l ) = {0,x- 1}. 

(3) . If s > 3 is odd then RJ w ({s} 1 ) = {0,x-l}. 
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(4) . Ifs is even then RJ s {{s} 1 ) = {0, i + (a; - l)/2, x - 1 : < i < I - 1}. 

(5) . If s> 3 is odd t/iera RJ(s, 1) = {0, a; — 1, a;}. 

(6) . J/s + t is odd i/ien RJ(s,t) = {0}. 

(7) . If s ^ t, s + t is even and t ^ 1 i/ien RJ(s, t) = {0, a; — 1}. 

(8) . If s ^= t, s + t is even and t ^ 1 i/ien RJ(s, t, s, t) — {0, a; — 1}. 

(9) . Ifr + s + t>5 is odd, r,s,t> 2, and r ^ t, then RJ(r,s,t) = {0}. 

(10) . Ifs is odd then RJ(r,s,r) = {0,x- 1}. 

(11) . If r + s + t > 6 is eijen, (r, s, i) ^ (4, 3, 5), (5, 3, A), and r,s, t ; > 2 are not aZZ i/ie 

same, then RJ(r, s, t) = {0}. 

(12) . i?Ji (2,l,l) = i?J(l,l,2) = i?J(4,3,5) = i?J(5,3,4) = {0,x- 1}. 

(13) . Ifs,t>0 and s + t is even then RJ w (l s ,2, 1*) = {0,x- 1}. 

(14) . //s > 1 i/ien .RJ'ioCl 3 , 2, I s , 2, I s ) = {0,35-1}. 

(15) . If s>2 is even then RJ w (l s , 2, l 8 " 1 , 2, l s+1 ) = {0,x- 1}. 

(16) . If s>2 is even then RJ w (l s+1 , 2, I s ' 1 , 2, I s ) = {0,x- 1}. 

(17) . J/ I is even then RJ{l l ) = {0, x - 1, x, 2x - 1}. 

(18) . RJ W {1, 2, 1) = {0, x - 1, 2.t - 1}. 

Remark 7.3. We are sure that we can remove the subscript 10 in RJ\q when more powerful com- 
putational tools are available to us. 

Proof. Even without the restriction of the bound 5, the inclusions RJ(s;p) C J(s\p) follow from 
Eq. (HI, Theorem PI [21, Theorems 2.14, 3.1, 3.5, 3.16, 3.18], and Theorem except in case 
(|7.2| and the last two cases. 

For s = 1', I > 2, Theorem [TT1 implies p- 1 S J(l i |p p( ^ 1) ) for all I > 1 and p > I + 3. In fact, 
by [HI Theorems 2.8, 2.14] we have 




(30) 



So if p > I + 3 then 





(mod p) if 2 { ?. 



(31) 



Further, setting ft., = H(l l ;p — 1), /io = 1 and h 



i = we have 



H(l l ;2p-l)=J2l 



E 



i 



i=Q l<ki<-<ki<p<k i + 1 <-<ki<2p 



kik 2 . . - h 




Here we have used geometric series expansion inside Q p such that for any m < p 




in 



h. 



-pX)ff((l i - 1 ,2,l wl - < );p-l) (modp 3 ) 



i=i 



= h m +p((m + l)h m+ i - hxh m ) 



h. 



+ (-l) m pH(m + l;p-l) (mod p 3 ) 
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by Eq. ((28]). It follows from equations (28]), (32]) and [ID Theorem 3.1] that 

1 z_1 

H(l l ; 2p-l) = 2hi + -J2 hih-i-i - (-l) l H(l;p- 1) (mod 



r i=0 

When I is even we have 

, , , (l + 2)H(l:p-l) 2H(l-l;p-l) , , 9n 

H(l l ;2p - 1) = 1 + y p[l _ x) 1 = ~2pB p + x (mod p 2 ). 

So it's divisible by p. This shows that 2x — 1 belongs to the reserved set RJ(l l ) in the penultimate 
case of the theorem. When I = 2n + 1 is odd hi = (mod p 2 ) and we get 

#(r; 2p - 1) = - } h 2l h 2n -2 l = = r~ (mod p) 

which is rarely congruent to 0. This explains why in case (2) we can't have 2x — 1 in the reserved 
set RJ(l l ) when I is odd. 

Finally let's turn to the last case of reserved set RJ(1, 2, 1). We have for any prime p > 7 

ff(l,2,l;2p-l) = V -L-=A + B + C + D, 



lm 2 n 

l<Km<n<2p 



where 



A= T T ^-=H(l,2,l;p-l) + -H(l,2;p-l), 

l<l<m<n<p 



D 



l<l<m<p<n<2p x r / 1 < fc<p f 

c = E rAr = f ff ( 1 ^ + E ? i ,1 i v 

* — ' lm z n \ vl * — ' p + m p + n 

l<l<p<m<n<2p v ^ 7 l<m<ri<p v ^ ' v ^ ' 

* = E rV= E 



.m 2 n ^-^ (p + Z)(p + m) 2 (p + n) 

p<l<m<n<2p l<Km<n<p v ^ ; v ^ ; 

Observe that D = ff(l,2,l;p- 1) = (mod p) by [HI Corollary 3.6] and H(l;p- 1) = (mod p 2 ) 
by Wolstenholme's Theorem. By geometric series expansion we get 

\ - 1 _ \ - (l-2p/m)(l-p/n) 2 
> 7 t^t r= > 7. (mod p ). 

e - — ' (p + mr[p + n) * — ' m z n 

l<m<n<p v ' y ' \<rn<n<p 

Hence 

ff(l,2,l;2p-l) = - (H(l, 2; p-l)+ #(2,l;p-l)) 
p 

-2if(3,l;p- 1) -H(2,2;p- 1) (mod p). 

By Theorem O and [2D Theorem 3.1] we have H(2, 2;p- 1) = #(3, l;p- 1) = (mod p). Further, 
from shuffle relation we have 

ff(l,2;p-l)+iy(2,l;p-l) = H (l;p - l)H(2;p - 1) - if(3;p - 1) = (mod p 2 ) 

by Theorem O This shows that H(l, 2, 1; 2p - 1) = (mod p). 

To prove the theorem we now only need to demonstrate that RJ(s; p) = J(s|p) for some p > |s| +3 
which can be done through a case by case computation. We put this part of verification online [2"U] . 
In fact, much more data are available in this supplement. □ 
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Are there any other s e N 1 (I < 3) besides those listed in Theorem 17.21 such that RJ(s) ^ {0}? 
In view of the last conjecture of [21] we believe there are. 

From Theorem 17.21 we see that to determine RJ(s) we often only need to study RJi(s) because 
for all non-homogeneous s not equal to (2r — 1. 1) or (1, 2, 1), the proportion of primes p such that 
Jt(s\p) — for all t > 2 is supposed to be positive. This implies that RJ(s) = RJi(s) for all such s. 
Precisely, we have the following 

Conjecture 7.4. Suppose s 6 N 1 such that (i) s = 1, or (ii) s = (1,2, 1), or (Hi) s = (2r — 1, 1) 
for some r > 1, or (iv) s = I 21 for some I > 1. Then RJ(s) — RJ2(s). For all other s we have 
RJ(b) = RJ 1 (b). 

When a > 2 we assume that H (s'; ap— 1) has random distribution modulo p 2 (the case a = 2 and 
s' = (1, 2, 1) has to be dealt with separately, but that's not hard). Thus the chance that p\H(s; ap) is 
less than 1/p 3 for large p. This implies that the probability of J2(s|p) = is roughly (1 — l/p 3 ) p — > 1 
as p — > oo. 

Definition 7.5. Let seN 1 . Define the reserved density of the MHS H(s) by 

density^s): X) = to^c P ■ N + 2 < p < X, J(s| P ) ^J(s)} 

Jy y ' tKprime p : |s| + 2 < p < X} v ; 

and the mth reserved density by 

density(i?J™(s); A) = « primfi P : N + 1 < P < ^> = RJ ™^ , (34 ) 

^ " ' j{prime p : |s| + 2 < p < X] v 7 



Conjecture 7.6. Let seN 1 . 77ien 

density (i? J (s); oo) 




i/Z = l,s>2, 

if 1 = 8 = 1 orl>2. 



Note that we always have RJ(s;p) C J(s\p). We have put the data strongly supporting Conjec- 
ture EH in 20, Appendix II]. In fact, we have only computed the first or the second reserved density 
because according to Conjecture 17.41 this is enough to determine the reserved density in whole. 

We now provide a heuristic argument for Conjecture 17.61 Suppose I = 1 and s = s > 2 first. 
Then by Proposition 16.11 we only consider H(s;r) for 1 < r < (p — 5)/2 + p(s — 1) because for 
most p the midpoint (p — l)/2 € J\(s\p) if and only if s is even (see Proposition 16. 3p . If we assume 
that when r varies the numbers H(s;r) distribute randomly modulo p for any large fixed prime p 
then the probability that J?(s\p) = RJi(s;p) is (1 - i/ p )(p-5)/2+p(s-i) _> as p oo. By 

Conjecture 17.41 (we also have a heuristic argument for it in this case) we see that the probability that 
J(s\p) = RJ(s;p) is lj\fe. 

Remark 7.7. Although we cannot prove the random distribution of H(s; r) for 1 < r < (p— 5)/2 + 
p(s — 1) modulo large prime p , in a recent paperj9], Garaev et al. show that for any s > 0, the set 
{H(s; r)) : r < p 1 / 2+£ } is uniformly distributed modulo a sufficiently large p. 

Now we assume I > 2. In general J°(s|p) does not have any symmetry so we see that the 
probability that J?(s\p) = i?Ji(s;p) is (1-1 /p) p - 5 -> 1/e asp — * oo, where 6 = t)i?Ji(s). When s docs 
not belong to the cases (i)-(iv) in Conjecture 17.41 we see that the probability that J(s\p) = RJ(s;p) 
is 1/e by Conjecture 17.41 

Finally let's deal with larger reserved sets when I > 2. By Theorem 17.21 we know that if s = l 2m 
ors= (1,2, 1) or s = (2r - 1, 1) for some r > 1 then RJ(s) = RJ 2 (s). Let s = l 2m . Then for 
p < n < 2p — 1 by definition H{\ 2m ;n) is a sum of many product terms with u p -value equal to either 
or 1, i.e., the denominator has at most one p factor. Assuming random distribution of pH (1 m ; n) 
modulo p 2 we see that when p < n < 2p — 1 the probability that p divides H(l 2m ; n) is 1/p 2 . Thus 
the probability that none of H(l 2m ; n) (p < n < 2p — 1) is multiple of p is equal to (1 — l/p 2 ) p ~ 2 . 
Similar heuristic argument implies that the probability that none of H(l 2m ;n) (2p < n < 3p) is 
multiple of p is equal to (1 — l/p 3 ) p , and so on. When Ip < n < (I + l)p and 2m < I < p the 
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probability is equal to (1 — l/p 2m+1 ) p . In general for p l < n < p t+1 (t > 1) we can break it into p 
subintervals and estimate inside each of the subintervals. It is easy to conclude that the probability 
that none of H(l 2m ; n) (p t < n < p t+1 ) is multiple of p is at least (1 - i/ p 2m+2t-iy*+\ Now 

^ 1 1 

p 2m+t-2 _ p 2m-2(p _ 1) ^ &S P ~~ * °°' 

We see that the probability that J(l 2m \p) = RJ(l 2rn ;p) is 

(1 - l/p) p ~ 2 — ► - asp^oo. 

e 

We omit the arguments for s = (1, 2, 1) and (2r — 1, 1) which are similar. 

We conclude our paper by some conjectures which concern distributions of irregular primes in 
disguised forms. 

Conjecture 7.8. Let r, s and t be positive integers. Then 

(1) . If s > 1 is odd then Ji(s\p) = {(p — l)/2,p — 1} for infinitely many primes p. 

(2) . If s > 1 is odd then Ji(s,s\p) = {(p — l)/2,p — 1} for infinitely many primes p. 

(3) . Let s,t be positive integers. Suppose s + t is odd. Then Ji(s,t\p) — {p — 1} for 

infinitely many primes p. 

(4) . Let r, s, t be positive integers such that r+s+t is odd and r =/= t. Then Ji(r, s, t\p) = {p — 1} 

for infinitely many primes p. 

Note that by various results of this paper and [21j an affirmative answer to any part of our 
Conjecture 17.81 would imply that there arc infinitely many irregular pairs (p,p — w) for any fixed odd 
number w (> 5 in case (|7.8jl ). Therefore, even if the sets of primes in Conjecture 17.81 are expected 
to be infinite they are extremely sparse; very likely they have zero density. 
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